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ABSTRACT 



We have identified and classified Zeeman pairs in the survey by Argon, Reid, &z 
Menten of massive star-forming regions with 18 cm ( 2 n 3 / 2 ,>/ = 3/2) OH maser emis- 
sion. We have found a total of more than 100 Zeeman pairs in more than 50 massive 
star-forming regions. The magnetic field deduced from the Zeeman splitting has allowed 
us to assign an overall line-of-sight magnetic field direction to many of the massive 
star-forming regions. Combining these data with other data sets obtained from OH 
Zeeman splitting, we have looked for correlations of magnetic field directions between 
star-forming regions scattered throughout the Galaxy. Our data do not support a uni- 
form, Galactic-scale field direction, nor do we find any strong evidence of magnetic field 
correlations within spiral arms. However, our data suggest that in the Solar neighbor- 
hood the magnetic field outside the Solar circle is oriented clockwise as viewed from the 
North Galactic Pole, while inside the Solar circle it is oriented counterclockwise. This 
pattern, including the magnetic field reversal near the Sun, is in agreement with results 
obtained from pulsar rotation measures. 



Subject headings: masers — ISM: magnetic fields — Galaxy: disk — H II regions - 
radio lines: ISM 
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1. Introduction 

Hydroxyl (OH) masers are an excellent tool for probing the environment of massive star- 
forming regions (SFRs). They are bright and easily observable and act as point tracers of the local 
velocity and magnetic fields, owing to their large Zeeman splitting coefficient. For the ground- 
state 2 n 3 / 2 , J = 3/2 transitions, the Zeeman splitting coefficient is as large as 0.59 km s _1 mG _1 . 
For a typical maser linewidth of 0.3 km s _1 , a milligauss-strength magnetic field is sufficient to 
cleanly separate Zeeman components of the 1665 MHz transition, and even fields as weak as several 
hundred microgauss can be detected. In addition, the sense of the splitting indicates the projected 
line-of-sight magnetic field orientation. 

Davies (1974) first hypothesized that a connection exists between magnetic field directions 
obtained from molecular Zeeman splitting and a large-scale Galactic magnetic field. He found 
that eight OH maser sources (and several H I clouds) mostly near the Solar circle with observable 
Zeeman splitting displayed magnetic fields aligned in the clockwise direction when viewed from the 
North Galactic Pole. Reid & Silverstein (1990) conducted a literature search and found a total 
of 17 reliable OH maser sources with detectable Zeeman pairs in the 2 n 3/ / 2 ,</ = 3/2 transitions. 
Of these 17, 14 had line-of-sight field directions consistent with an overall clockwise Galactic field. 
In the largest survey to date, Baudry et al. (1997) extended this dataset to include a total of 46 
sources with Zeeman pairs in the 2 n 3 / 2 , J = 3/2 and J = 5/2 transitions of OH. Excluding the 
Galactic center, they found 28 sources consistent with a clockwise magnetic field and 17 sources 
consistent with a counterclockwise field. 

Argon, Reid, & Menten (2000, henceforth Paper I) conducted a survey of 91 interstellar 
regions with known 2 n 3/2 , J = 3/2 OH maser emission. They identified and mapped maser spots 
in the regions. In this paper we extend the analysis of Paper I by identifying Zeeman pairs toward 
massive SFRs, many of which contain ultracompact (UC) H II regions. We also investigate whether 
or not the magnetic fields deduced from OH Zeeman splitting around massive SFRs are related to 
the large-scale magnetic field of the Galaxy. 

2. Observations 

A total of 91 maser sources were observed by Paper I with the A configuration of the National 
Radio Astronomy Observatory's Very Large Array 1 in Socorro, NM. All sources were observed 
in both left- and right-circular polarization in at least one main- line 2 n 3/2 , J = 3/2 18 cm OH 
transition (1665.4018 and 1667.3590 MHz), and some were also observed in one or both of the 
satellite-line transitions (1612.2310 and 1720.5300 MHz). The 0.763 kHz spectral channel separation 
used corresponds to a velocity separation of 0.14 km s _1 for the main-line transitions (FWHM 



lr The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under 
cooperative agreement by Associated Universities, Inc. 
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~ 20% larger). Further details of the observations can be found in Paper I. 

3. Zeeman Pair Selection Criteria 

In the presence of a magnetic field, the degeneracy of the magnetic sublevels of the OH molecule 
is broken. For any given hyperfine transition, this results in three or more spectral lines with 
differing polarization properties — the Zeeman effect (see Davies 1974, for more details). OH 
masers often display pairs of strong left and right elliptically polarized a-components. For reasons 
that are not well understood, purely linearly polarized 7r-components are rarely, if ever, detected. 

The task of identifiying Zeeman pairs — two oppositely elliptically polarized lines at the same 
position on the sky — is relatively simple, provided one has sufficient angular resolution to resolve 
maser clusters (e.g., Reid et al. 1980). However, with the VLA one can only achieve an angular 
resolution of « 1 arcsec for the OH transitions at wavelengths of 18 cm. Were the emission in any 
spectral channel to come from only one (unresolved) spot on the sky, this would not be a significant 
limitation. However, in reality the emission in any given spectral channel can come from spectral 
blends of masers at different positions in the source separated by less than the beam size. Thus, 
great care must be taken to determine Zeeman pairs from the Paper I data. 

In order to identify Zeeman pairs, a modified version of the NRAO Astronomical Image Pro- 
cessing System (AIPS) task JMFIT was used to identify the brightest maser emission spot in each 
spectral channel (velocity plane) for any given source. A maser line was identified when emission 
was seen from approximately the same position in at least three adjacent velocity channels, with 
the strongest emission not found on the edge channels. This resulted in a large data set that could 
be automatically searched for Zeeman pairs. 

Ideally, one would like to find Zeeman pairs in which the oppositely circularly polarized com- 
ponents were coincident on the sky to better than their apparent spot sizes. However, Reid et al. 
(1980) showed that OH masers in W3(OH) clustered on a scale of 3 x 10 15 cm (0"09 at 2.2 kpc; 
distances are discussed in §4.1 and listed in Table 1) and that oppositely circularly polarized maser 
lines within a cluster gave reliable magnetic field measurements when interpreted as Zeeman pairs. 
Thus, we searched for oppositely circularly polarized maser lines with projected positions separated 
by less than 3 x 10 15 cm. While the main-line 1.6 GHz cr-components are in general elliptically 
polarized, it is almost functionally equivalent to treat them as being circularly polarized for the 
purposes of right-circular polarization (RCP) and left-circular polarization (LCP) VLA data. The 
caveat is that, for example, a right-elliptically polarized maser component, being the sum of circu- 
larly and linearly polarized emission, may appear in both RCP and LCP data at the same velocity, 
although the LCP emission will be much weaker, as the LCP feed is sensitive to only half of the 
(presumed small) linear component of the full elliptical polarization. 

Some Zeeman pairs consisted of single RCP and LCP lines that were coincident on the sky 
to better than 3 x 10 15 cm. We consider these easily recognizable and reliable Zeeman pairs. 
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However, more than one maser line in the same polarization would often be seen within the assumed 
3 x 10 15 cm clustering length, such that it was unclear which should be paired with the maser spot 
in the other polarization. While such data are insufficient to make an unambiguous Zeeman pair 
identification in all such cases, sometimes corroborating evidence from another maser transition 
(or occasionally even the same transition) in the same spatial region allowed the identification of 
Zeeman pairs. In other cases, the line-of-sight field direction was deducible, even if unambiguous 
identification of a Zeeman pair was not possible, for example, when multiple LCP and RCP lines 
were seen, but each RCP line was at a higher velocity than any of the LCP lines. We felt it 
important to generate specific, objective, criteria that allowed an algorithmic approach to search 
for Zeeman pairs and accept or reject them without resorting to a visual examination of spectra 
on a case-by-case basis, which can lead to significant unquantifiable biases. Below we describe how 
we handled each of the complications. 

For any given source, a "window" of radius 3 x 10 15 cm was centered on the strongest line 
in either polarization in a single transition. In principle, if the window included one line in each 
polarization, the lines would be identified as a Zeeman pair. However, emission channels key to 
identifying the second maser line in a Zeeman pair sometimes lay just outside the window when it 
was centered on the position of brightest emission, whereas were the window moved between the 
position of the lines, it would pass the < 3 x 10 15 cm separation criterion. To overcome this problem, 
the window center was moved from the centroid of strongest emission along a grid, in increments 
of 5% of the window radius, out to a maximum displacement of 50% of the radius in both the R.A. 
and Dec directions. The window containing the largest number of lines was used, provided that at 
least one line of each polarization was included. "Tiebreaker" criteria consisted first of choosing 
the window with the largest number of velocity channels represented and second of choosing the 
window whose offset from the grid center was smallest. Next, all velocity channels contained in the 
chosen window were flagged as used, and the process was iterated with the remaining maser spots. 
In this way, the entire source was partitioned into windows consisting of potential Zeeman pairs. 
This partitioning was done independently for each maser transition for which data were collected 
for the source. 

Identification of Zeeman pairs in windows containing only one maser line in each polarization is 
straightforward. Figure la shows spectra from two windows for G351. 775— 0.538, one in each of the 
1665 and 1667 MHz transitions. The 1665 MHz window contains two oppositely circularly-polarized 
maser lines. Removing the effects of the Zeeman splitting from a magnetic field of —3.80 ±0.03 mG 
(i.e., oriented in the hemisphere toward the Sun.), the center velocity of the emission would be at 
—26.74 ± 0.02 km s _1 . The 1667 MHz window, centered 120 mas away, shows a similar splitting 
of —3.66 ± 0.06 mG centered at a velocity of —26.79 ± 0.02 km s _1 . While a large fraction of 
the Zeeman pairs we identified were found in windows containing exactly one maser line in each 
polarization, a few of these pairs appeared with a corroborating nearby pair of a similar magnetic 
field in another transition. 

Table 1 is a list of Zeeman pairs with "non-zero" velocity separation. The flux, velocity, and 
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velocity width of lines was determined by a Gaussian fit to the center 5 channels of each line when 
available, else to 4 or only 3 channels. (By our criteria for a maser line, the center 3 channels must 
necessarily be detected.) Pairs are only included in this table if the velocity difference between the 
LCP and RCP lines exceeds twice the signal-to-noise-limited velocity error, taken to be the sum 
in quadrature of the HWHM line width divided by the signal-to- noise ratio of each of the lines. 
However, in no case is a pair of lines allowed to satisfy this criterion if the velocity center separation 
is less than 5% of the quadrature sum of the line widths, or typically about 0.02 km s _1 . These 
"pairs" may be actual Zeeman pairs with a small splitting due to a small magnetic field, or they 
may be single maser spots with a large linear polarization fraction, but the limited spatial and 
velocity resolution of the VLA and the absence of full-polarization observations render it difficult 
to decide between these two cases. These "near zero" velocity separation maser spots are listed in 
Table 2. 

While identification of many of the pairs listed in Table 1 was straightforward, not all windows 
were so simple to interpret. When a window contains more than one maser line in a polarization, 
interpreting the Zeeman pattern is more complicated. In this case, it may still be possible to 
know the light-of-sight direction of the magnetic field, but not its strength. Figure lb shows the 
spectrum of one window associated with G032. 744— 0.076. There is one RCP maser line (labelled 1) 
centered at 32.38 km s _1 and two LCP lines (2 and 3) at higher velocities (34.60 and 35.02 km s _1 , 
respectively). Additionally, there is one feature (4) that does not meet our criteria as a maser line, 
since the low-velocity side of the line is not seen, probably owing to blending with a stronger line 
elsewhere in the source. It is unclear which LCP line should be associated with the RCP line as 
a Zeeman pair. Regardless of the exact pairing, the RCP line is at lower velocity than any LCP 
line, so the line-of-sight orientation of the magnetic field is known to be pointed toward the Sun. 
Cases such as this, where the field direction is known but its strength is not, are frequent among 
our data. When a range of values is quoted for the magnetic field in Table 1, the direction of the 
magnetic field is well-determined, and the smallest and largest possible magnetic field intensities 
are given. Some line parameters not listed in the table are given in the notes on individual sources. 

Sometimes when a window contains enough maser lines such that the magnetic field strength 
or direction is ambiguous, supporting evidence may exist to suggest one pairing above all others. 
If a pairing in one transition implies a magnetic field strength and direction consistent with that 
of a pair in another transition at nearly the same location, that pairing is accepted in preference 
to all other possible pairings. Figure lc shows a 1665 MHz and a 1667 MHz window in K3— 50, 
coincident to within 46 mas. In each window, there is one RCP line (1) at lower velocity than 
either of the two LCP lines (2 and 3), so the line-of-sight field is oriented toward the Sun. But 
the field magnitude can also be deduced, assuming that the nearly overlapping windows in the two 
transitions are measuring the field at approximately the same location. In the 1665 MHz window, 
pairing lines 1 and 2 would give a magnetic field of —1.5 mG, while the 1-3 pairing implies a —2.5 mG 
field. In the 1667 MHz window, the 1-2 pairing yields —0.8 mG, while 1-3 yields —2.6 mG. The 
close agreement between the field values for the 1-3 pairs for the two transitions suggests that 
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these Zeeman pairs should be favored over the other pairing possibilities. For definiteness, when 
appealing to corroborating evidence from one transition to resolve pairing ambiguities in another 
transition, we required AB < 0.4 mG, Av c < 0.2 km s~ 1 , and A9 < 0"2sec2:tanz, where AB is 
the difference of the implied magnetic field strengths, Av c is the separation of the center velocity 
of each pair, AO is the positional offset between pairs, and z is the zenith angle of the source as 
viewed from the location of the VLA at transit. As noted in Paper I, relative positional offsets 
between different OH transition frequencies are ~ O'fl for high-declination sources, but grow to 
several arcseconds for sources at very low declination (< —40°). The function "sec z tan z" was 
chosen in the spirit of the plane-parallel approximation for the atmospheric contribution to the 
relative error between observations of different transitions (see Reid et al. 1999). 

If two pairings in the same window in the same transition imply the same magnetic field 
direction and strength, those pairings are also accepted over all other pairings. This situation 
arises when two Zeeman pairs in the same transition are encompassed in the same window. Figure 
Id shows a 1665 MHz window for G45. 465+0. 047 that contains a cluster of maser lines. Since all 
of the RCP lines are at a higher velocity than the LCP lines, the magnetic field is oriented away 
from the Sun. However, the velocities associated with the line centers favor two particular pairings: 
pairing lines 1 and 4 yields a magnetic field strength of 3.7 mG, and 2 and 5 yield 3.8 mG. 

Occasionally, a window would produce a spectrum with a clear, strong line in one polarization 
and an incomplete line in the opposite polarization. A closer look at the data would show that 
one or two critical points of the incomplete line lay just outside the window. By increasing the 
size of the window, a Zeeman pair was often identified. Using a larger window (generally twice 
the clustering radius) can be justified for two reasons. First, the window radius of 3 x 10 15 cm 
is based on the clustering scale analysis of Reid et al. (1980), which is itself based on only one 
source: W3(OH), assuming the photometric distance of 2.2 kpc (Humphreys 1978). The kinematic 
distance of 4.3 kpc, assuming the NH3 LSR velocity of —45 km s _1 (Reid, Myers, & Bieging 1987) 
and a flat Rq = 8.5 kpc, ®(R) = 220 km s _1 rotation curve, would correspond to a clustering 
scale of 6 x 10 15 cm, or twice that which is derived from the photometric distance. In any case, 
it is unclear to what extent this varies among massive SFRs. Second, the distance to some of the 
massive SFRs is unknown. If the assumed distance to the source is greater than the actual distance, 
the (angular) window radius used would be proportionally smaller than the desired 3 x 10 15 cm. 
For sources for which the kinematic distance ambiguity has not been resolved, the Zeeman pair 
identification analysis was done twice, assuming each of the near and far distances for calculation of 
the clustering scale. Zeeman pairs identified assuming either distance were accepted and are listed 
in Table 1. Despite the distance uncertainty, we believe that Zeeman pairs identified around these 
sources are still reliable indicators of the local magnetic field. When the far kinematic distance 
is more than a few times the near kinematic distance, possible errors involve partitioning a far 
source into windows that are too large or a near source into windows that are too small. In the 
former case, nearly all of the maser emission of a source (presumedly) located at the far kinematic 
distance would be concentrated in a single window when the near distance was used for the angular 
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clustering scale, resulting in confusion of Zeeman pairs. In the latter case, a (presumed) near source 
would be micropartitioned into windows of small radius, rejecting many Zeeman pairs coincident 
to within the desired linear clustering scale. Thus, misresolving the near/far kinematic distance 
ambiguity results in the identification of fewer Zeeman pairs, regardless of whether a far source is 
assumed to be near or a near source is assumed to be far. 

Zeeman pairs were also graded to reflect the uncertainty in pairing. Pairs from windows 
containing only one line in each polarization were assigned a grade of B, unless they were confirmed 
by a similar pair in a second transition, in which case they were given a grade of A. All pairs 
whose pairing was ambiguous were given a grade of C, even if the ambiguity was resolved by 
appealing to a second pair in the same or another transition. Pairs for which the separation of the 
Zeeman components exceeded the default window size were assessed a grade of D. In one case in 
G081. 871+0. 781, a grade of D was assigned to a Zeeman pair believed to be spurious; for further 
details see §4.3. This grading system is in the spirit of, although not identical to, the system used 
in Reid & Silverstein (1990), in which the main emphasis was on whether interferometry was used 
to identify Zeeman splitting, and if so, at what resolution. All of the Zeeman pairs listed in Table 
1 were identified using interferometry. 
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4. Zeeman Pairs 
4.1. Distances 

The distances given in Tables 1 and 2 are kinematic distances assuming a flat rotation curve 
with Rq = 8.5 kpc and Q(R) = 220 km s -1 , except as noted in the text. For sources inside the solar 
circle in the first and fourth Galactic quadrants, there is a distance ambiguity. For a number of these 
sources, the kinematic distance ambiguity is resolved from H I absorption observations by Fish et al. 
(2003). That paper also provides an analysis of the scale height of known UCH II regions (embedded 
in massive SFRs) above the Galactic plane, including an analysis of the accuracy of resolving the 
distance ambiguity to a UCH II region given only the kinematic distances and Galactic latitude of 
the region. When the distance was not known by any other method, the Galactic latitude was used 
to resolve the kinematic distance ambiguity toward massive SFRs only when the difference between 
the near and far kinematic distances exceeds 4.5 kpc, corresponding to 90% predictive accuracy 
(Fish et al. 2003). The majority of these sources are located in the fourth Galactic quadrant. 

When evaluating kinematic distances, we assumed that massive SFRs have a peculiar motion 
of < 10 km s _1 . This freedom is important to placing a few sources (e.g., G355. 345+0. 146) at a 
sensible distance. Nevertheless, it is possible for peculiar motions to exceed 10 km s -1 . Hydroxyl 
and H2CO absorption measurements toward G10. 624— 0.385 and other sources in the W31 complex 
demonstrate deviations from the velocities expected from the rotation model of as much as 36 
km s~ l (Wilson 1974). 

For sources whose median LSR velocity of OH maser emission is close to km s _1 , the near 
kinematic distance is often unrealistically close. Since the angular window radius as described in 
§3 is inversely proportional to the assumed distance to the source, a source assumed to be much 
nearer than it is would be partitioned into windows encompassing too large of an angular area. The 
resulting windows would be few in number and contain a large number of maser lines. The effect 
is amplified for especially near distances (a few hundred parsecs or less), where a small velocity 
error translates into a large fractional distance error and therefore a large fractional error in the 
window size. To be consistent (if somewhat arbitrary) in handling these extremely near sources 
(<S 2 kpc), we used the distance produced by averaging the kinematic distance derived from the 
median OH maser velocity and 2.0 kpc. Since massive SFRs are few and far between, it is unlikely 
that a massive SFR whose distance is not independently known is actually located within 1 kpc of 
the Sun. 



4.2. Satellite-Line Field Splitting 

Unlike the main-line (1665 and 1667 MHz) transitions, which split into only one pair of u- 
components, the satellite-line transtitions split into three pairs with line intensities in the ratio 
6:3:1 (Davies 1974). As mentioned in a footnote to Table 1, the Zeeman splitting coefficient 
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8.2 


1665 


4.66 


-51.83 


0.48 


2.60 


-52.11 


0.63 


0.032 


-0.5 


B 



a Kincmatic distances assuming i?o = 8.5 kpc and 0o = 220km s 1 , unless note in text (§4.3) says otherwise. 

b Pcak flux density, center velocity, and line width (FWHM) of left- and right-circular polarized lines, respectively, as determined by Gaussian fit. 

c Angular offset between LCP and RCP lines. When a range of values is given for the magnetic field (B), the offset is stated for the lines of peak flux density in 
each polarization. 

d Full magnitude of B field. The sign indicates the direction of the line-of-sight component (positive pointing away from Sun). The ratio of the velocity separation 
between RCP and LCP components to B field was taken to be 0.590, 0.354, 0.236, and 0.236 km s _1 mG -1 for the 1665, 1667, 1612, and 1720 MHz transitions, 
respectively. See §4.2 for further details. When a range is given, the extrema represent the minimum and maximum field strength implied by all possible pairings 
of RCP and LCP lines. 

e See text (§3) for description. 

'Likely spurious. See §4.3 for details. 
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for the satellite-line (1612 and 1720 MHz) transitions was assumed to be 0.236 km s _1 mG -1 . 
This splitting coefficient is correct for the intensity-weighted average of the three lines in each 
polarization. However, maser amplification may be greatest for the strongest line in each triad. If 
only the strongest line is seen after amplification, the splitting coefficient would be 0.114 km s _1 
mG" 1 , or approximately half of the intensity- weighted value. Indeed, the splitting coefficient could 
be between the two extremal values, as would be consistent with amplification of all three o- 
components in a partially-saturated maser. A comparison of the magnetic field strengths inferred 
from 1720 MHz pairs in W3(OH) listed in Table 1 with those shown on VLBI maps (Bloemhof, Reid, 
& Moran 1992) is more consistent with a splitting coefficient closer to 0.236 km s _1 mG -1 . Analysis 
of field strengths inferred from 1612 and 1720 MHz Zeeman pairs in all sources also supports a 
large splitting coefficient (see §4.4). In any case, the direction of the line-of-sight component of the 
field is unaffected. 



4.3. Notes on Specific Sources 

G000. 670— 0.058: The separation between the LCP and RCP components of the pair exceeds 
the assumed 3 x 10 15 cm correlation length by less than a factor of two. 

G005. 886— 0.393: The two pairs receiving an A grade are separated by approximately 0".3. 

G009. 622+0. 195: We adopt a distance of 5.7 kpc to this source, consistent with observations 
by Scoville et al. (1987) and Hofner et al. (1994), although the distance obtained kinematically 
could be as near as 0.7 kpc. The second Zeeman pair listed in Table 1, offset by about 5" from the 
other two pairs, is in a separate maser site within the same massive star-forming cluster. Additional 
lines for the first pair in table 1 receiving a grade of C consist of an LCP line of 1.65 Jy at —3.97 
km s _1 and an RCP line at 7.85 Jy and 1.22 km s _1 . For the second such pair, there are four 
additional LCP lines: 0.79 Jy at 7.04 km s _1 , 0.39 Jy at 6.50 km s _1 , 0.47 Jy at 6.09 km s~\ and 
1.32 Jy at 5.49 km s -1 . 

G012. 680— 0.181: Additional lines for the 1665 pair receiving a grade of C consist of two RCP 
lines: 2.23 Jy at 60.71 km s _1 and 2.58 Jy at 61.17 km s _1 . For the 1667 pair, there are two 
additional LCP lines of 2.92 Jy at 64.81 km s" 1 and 2.28 Jy at 62.76 km s" 1 and one additional 
RCP line of 3.30 Jy at 62.15 km s _1 . From observations using the 43 m telescope at Green Bank, 
Zheng (1991) reports two Zeeman pairs of around —5 mG each. The 4 kpc distance assumed for 
W33 B and W33 A is the kinematic distance taken from Haschick & Ho (1983). 

G020.081-0.135: There is an additional LCP line of 4.44 Jy at 45.47 km s" 1 for the pair 
receiving a C grade. 

G032. 744-0. 076: There is an additional RCP line of 0.24 Jy at 40.24 km s" 1 for the pair 
receiving a C grade. 

G034.257+0.154: There is an additional LCP line of 2.24 Jy at 57.97 km s" 1 for the pair 
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receiving a C grade. This source is one of three in the G34.3+0.2 complex. The magnetic field 
structure of this cometary H II region, source C in Zheng, Reid, & Moran (2000), is unclear. The 
magnetic field of source B (G034. 258+0. 153) is ordered and pointing primarily away from the Sun. 

G035.024+0.350: There is an additional LCP line of 0.54 Jy at 43.83 km s" 1 and an RCP line 
at 1.83 Jy and 44.63 km s _1 . 

G035.197-0.743: There is an additional LCP line of 1.04 Jy at 26.35 km s _1 . 

G035. 577— 0.029: In each of the two pairs, the angular offset between the LCP and RCP lines 
is consistent with a separation of about 2 x 10 16 cm, assuming the kinematic distance. There are 
two additional LCP lines of 2.85 Jy at 51.89 km s _1 and 2.61 Jy at 49.91 km s _1 for the second 
pair receiving a D grade. 

G040.622-0.137: There is an additional RCP line of 1.04 Jy at 34.76 km s" 1 for the 1665 pair 
receiving a grade of C. For the 1667 pair, there are two additional LCP lines: 2.82 Jy at 31.53 
km s _1 and 1.34 Jy at 30.85 km s _1 . 

G045.071+0.134: For the 1665 pair, there are two additional LCP lines of 6.24 Jy at 53.57 
km s _1 and 8.53 Jy at 52.95 km s _1 and one additional RCP line of 7.21 Jy at 55.91 km s^ 1 . For 
the first 1667 pair, there is an additional RCP line of 7.21 Jy at 55.91 km s _1 . For the second 1667 
pair, there is an additional LCP line of 2.97 Jy at 55.78 km s _1 and an RCP line of 0.91 Jy at 
57.17 km s _1 . 

G045. 465+0. 047: The LCP and RCP lines of each of the two pairs are located in a circle of 
3 x 10 15 cm. While the pairing of these lines is theoretically ambiguous, we feel that the identical 
magnetic fields implied by the listed pairings and the small velocity difference between these pairings 
(0.65 km s _1 , less than the 0.8 km s _1 turbulent velocity differences expected within a masing cloud) 
are sufficient to justify our claimed magnetic field values. 

G049.469-0.370: In 18-cm VLA observations, Gaume & Mutel (1987) also find a Zeeman pair 
with a positive magnetic field. For more information, see the notes for the next source. 

G049.488-0.387: Gaume & Mutel (1987) find about half a dozen Zeeman pairs in the W51 
complex around regions el and e2, all indicating a positive magnetic field. In 5-cm VLBI observa- 
tions, Desmurs & Baudry (1998) find a Zeeman pair between regions el and e2. However, using 
the VLBA Argon, Reid, & Menten (2002) find two very strong 1665 MHz Zeeman pairs indicating 
magnetic fields of —21 and —19.8 mG. The kinematic distance is taken from Caswell (2001). For 
the 1665 pair, there is an additional LCP line of 59.42 Jy at 58.90 km s _1 . 

G069. 540— 0.976: For the 1665 pair receiving a C grade, there are two additional LCP lines: 
5.07 Jy at 15.45 km s _1 and 3.17 Jy at 15.06 km s _1 . Preliminary (as yet unpublished) VLBA 
18-cm observations by the authors of this paper and VLBI 5-cm maps of Desmurs & Baudry (1998) 
confirm that the magnetic field is consistently oriented toward the Sun across the source. 

G070. 293+1. 601: For the 1665 pair receiving a C grade, there is an additional LCP line of 
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4.03 Jy at -20.39 km s _1 . For the 1667 pair, there is an additional LCP line of 2.46 Jy at -20.80 
km s _1 . However, the two grade-C pairs in the table meet the self-consistency criteria (AO < 0".2, 
Av c < 0.2 km s _1 , AB < 0.4 mG between pairs in the two transitions). 

G081. 721+0. 571: Both this source and G081. 871+0.781, believed to be in the same complex, 
have average LSR velocities of a few km s" 1 positive. This puts these sources near the tangent 
point at this Galactic longitude with a large uncertainty in distance. We adopt the Dickel, Dickel, 
& Wilson (1978) distance of 2.0 kpc as the distance to the W75 complex. For the pair receiving a 
grade of C, there are additional LCP lines of 15.48 Jy at 0.93 km s _1 , 8.45 Jy at 0.58 km s _1 , and 
2.44 Jy at -0.27 km s" 1 as well as RCP lines of 0.40 Jy at 4.79 km s -1 , 0.65 Jy at 3.53 km s _1 , 
and 1.49 Jy at 2.97 km s _1 . 

G081. 871+0. 781: Preliminary (as yet unpublished) VLBA observations indicate a reversal of 
the line-of-sight magnetic field direction across the source. The —25.4 mG pair is not seen with 
the VLBA. However, a cluster that contains a large number of maser spots with a large velocity 
spread (pa 20 km s _1 ) is detected. It is probable that this "pair" is actually two nearby lines whose 
velocity difference is due to a velocity gradient rather than an extremely large magnetic field. For 
distance information, see the note for G081. 721+0. 571. 

G109. 871+2. 114: For the 1665 pair receiving a C grade, there is an additional LCP line of 
13.65 Jy at —13.90 km s -1 . However, the ambiguity is resolved by examining the 1667 transition, 
which is within 0".2 of the 1665 pair in question. Cohen, Brebner, & Potter (1990) find a Zeeman 
pair that implies a magnetic field direction oriented away from the Sun. Monitoring of this Zeeman 
pair over time shows that the strength of the magnetic field is decreasing. 

Gill. 533+0. 757: The distance of 2.8 kpc is taken from Campbell & Thompson (1984). 

G133.946+1.064: For the 1612 pair, there is an addition RCP line of 3.76 Jy at -42.92 km s -1 . 
For the 1720 line receiving a C grade, there is an additional LCP line of 6.88 Jy at —45.23 km s _1 
and two additional RCP lines: 4.47 Jy at —42.72 km s _1 and 7.47 Jy at —44.36 km s -1 . Previous 
VLBI observations have clearly established that the magnetic field is predominantly oriented away 
from the Sun (Garci'a-Barreto et al. 1988). The 2.2 kpc distance used for W3(OH) is based on OB 
star luminosities from Humphreys (1978). 

G343. 128-0.063: There is an additional LCP line of 54.52 Jy at -33.93 km s _1 . 

G344.581-0.022: The near kinematic distance of 0.6 kpc (Forster & Caswell 1989) is too 
close, as evidenced by the cluttered spectra produced in circles of radius 3 x 10 15 cm, but no pairs 
are seen at all at the far kinematic distance of 15.8 kpc. A distance of 2.0 kpc was used, which 
produces circles of radius 10 15 cm assuming the near kinematic distance is correct. For the first 
pair receiving a C grade, there is an additional LCP line of 1.99 Jy at —4.66 km s _1 and an RCP 
line of 2.93 Jy at —4.01 km s -1 . For the second such pair, there is an additional LCP line of 3.74 
Jy at —1.16 km s -1 . 

G345. 011+1. 792: For the pair receiving a C grade, there is an additional LCP line of 3.39 Jy 
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at -21.67 km s _1 . 

G345. 505+0. 347: For the pair receiving a C grade, there is an additional LCP line of 2.72 Jy 
at -16.92 km s -1 . 

G348. 549— 0.978: For the pair receiving a C grade, there are two additional RCP lines: 4.46 
Jy at -13.11 km s" 1 and 2.64 Jy at -19.26 km s _1 . 

G351. 416+0. 646: For the pair receiving a C grade, there is an additional LCP line of 22.29 Jy 
at -10.33 km s -1 . 

G351. 775— 0.538: For the 1665 pair receiving a C grade, there is an additional RCP line of 
3.04 Jy at —10.72 km s -1 . For the 1667 pair receiving a C grade, there is an additional RCP line 
of 1.11 Jy at 0.25 km s _1 . Preliminary VLBA observations indicate a reversal of the line-of-sight 
magnetic field direction across the source (Argon, Reid, & Menten 2002). 

G353.410-0.361: There is an additional LCP line of 3.55 Jy at -20.56 km s^ 1 and an RCP 
line of 2.34 Jy at —19.12 km s _1 . Caswell (2001) finds evidence for a magnetic field pointing in 
the opposite direction based on 6030- and 6035-MHz OH masers. 

G355. 345+0. 146: The rotation curve is quite flat in this part of the galaxy. The distance of 
23.1 kpc was derived assuming a velocity of +7.5 km s _1 , or 10 km s^ 1 less than the actual average 
value as deduced from OH maser emission. (Recall that a positive velocity in the fourth Galactic 
quadrant indicates that the source is located outside the solar circle, or at least 17 kpc away at 
this Galactic longitude.) The distance deduced from the actual average velocity of +17.5 km s" 1 
is unphysically large, assuming a flat rotation curve with Rq = 8.5 kpc and Q(R) = 220 km s -1 . 
There are two additional LCP lines: 15.38 Jy at 18.95 km s^ 1 and 15.84 Jy at 18.55 km s -1 . 

4.4. Distribution of Magnetic Field Strengths 

A histogram of the inferred magnetic field strengths for 75 Zeeman pairs from Table 1 is 
included in Figure 2. Only those Zeeman pairs implying a single value of the magnetic field strength 
are included, and the —25.4 mG pair in G081. 871+0. 781 is omitted (see §4.3). A Gaussian fit is 
superimposed on the histogram. The distribution of field strengths has a HWHM of 3.9 ±0.1 mG 
and is centered at mG to within the formal error. This is consistent with average field magnitudes 
of several milligauss deduced from OH Zeeman splitting in other studies (e.g., Garci'a-Barreto et 
al. 1988). Note that this field strength is what would be expected for magnetic field strength 
enhancement of \B\ oc n 1 / 2 (Mouschovias 1976) from interstellar values of B and n (~ 1 /j,G, 
~ 1 cm~ 3 ) ton~ 10 6 cm~ 3 , as is typical for sites of OH masing (Reid & Moran 1981). 

Of the 15 Zeeman pairs detected at 1612 and 1720 MHz, 10 imply a magnetic field strength 
less than or equal to the HWHM assuming a splitting coefficient of 0.236 km s _1 mG -1 , 3 imply a 
field strength greater than the HWHM, and the 2 pairs with a grade of C in Table 1 are inconclusive 
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owing to pairing ambiguities. If a splitting coefficient of 0.114 km s _1 mG" 1 is assumed, 3 Zeeman 
pairs imply a magnetic field strength less than the HWHM and 12 imply a greater field strength. 
Thus, magnetic field strengths deduced from 1612 and 1720 MHz Zeeman pairs are more consistent 
with values obtained from the main-line transitions when a splitting coefficient value closer to the 
high value (0.236 km s -1 mG' 1 ) is assumed. While we cannot rule out the possibility that the 
physical requirements for masing at 1612 and 1720 MHz cause satellite-line maser spots to be 
observed preferentially at locations of increased magnetic field strength, we feel justified in using 
the larger splitting coefficient. 



4.5. Other OH Zeeman Studies 

This study was undertaken primarily to investigate Davies's hypothesis that the magnetic 
fields deduced from OH Zeeman splitting in massive SFRs exhibit Galactic-scale organization. 
Preliminary analysis of the data from the Paper I survey yielded 40 massive SFRs for which OH 
Zeeman splitting suggested a local field direction (Fish et al. 2002). While the data set was 
homogeneous, it excluded southern sources, which are not observable from the VLA. To address 
this, the current study includes OH Zeeman pairs from a number of other studies as well. Chief 
among these are the 6 GHz surveys of Baudry et al. (1997) and Caswell &; Vaile (1995), the 
1.6 GHz study of Gaume & Mutel (1987), and the original literature search of Reid & Silverstein 
(1990). Two sources for which we have unpublished maps in the 1665 and 1667 MHz OH transitions 
with the VLBA have also been included (ON2 N and S269) (in preparation). Other supplementary 
papers relevant to only one or a small number of our sources are mentioned with the appropriate 
source(s) in §4.3. The objective is to compile a list of massive SFRs restricted to those with a 
predominant magnetic field direction. 

When assigning an overall magnetic field direction to a massive SFR, we require that at least 
two-thirds of all Zeeman pairs for any one source imply the same line-of-sight field direction. In 
addition to problems associated with a possible magnetic field reversal across the source, the VLA 
beam is large enough to encompass several maser clumps. While this could potentially blend maser 
spots, resulting in misidentification of Zeeman pairs, we have striven to obtain a sound collection 
of Zeeman pairs by discarding all pairs that do not meet the objective criteria outlined in §3, at 
the expense of completeness. The even lower spatial resolution of the data taken from A = 5 cm 
studies is not a concern since the intensity of RCP and LCP lines in a 6035 MHz pair is typically 
much more equal than is seen in ground-state Zeeman pairs (Caswell & Vaile 1995). This fact, 
combined with a much smaller Zeeman splitting coefficient, makes unambiguous pairing of maser 
spots easier in the 6 GHz transitions. Thus we do not believe that blending of multiple maser spots 
in the same beamwidth is a significant source of contamination for the magnetic field measurements 
in this study. The final sample to be used in §5.2 comprises 45 of the 53 sources in Table 1, plus 
29 from other information. 
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5. Discussion 

A number of historical findings motivated us to study whether or not the Galactic magnetic 
field is sensed by OH maser Zeeman splitting. The original hypothesis of a link between the 
Galactic magnetic field and interstellar OH masers was advanced by Davies (1974), who, noting 
that his compilation of eight OH maser sources all displayed a magnetic field consistent with the 
clockwise direction of Galactic rotation, suggested that the OH masers traced a circular Galactic 
magnetic field. For the most part, the OH sources included in his study are located near the 
solar circle, precluding investigation of reversals of the large-scale magnetic field direction with 
Galactocentric radius. Reid & Silverstein (1990) tested the Davies hypothesis by finding 12 new 
OH maser sources, of which 10 had a line-of-sight magnetic field aligned with the direction of 
Galactic rotation. This seemed to be consistent with Davies's suggestion and a concentric-ring 
model of the Galactic magnetic field as deduced from the pulsar rotation measure study of Rand &; 
Kulkarni (1989). However, Rand k, Kulkarni also presented evidence for magnetic field reversals 
between rings, which was not critically tested in the Reid & Silverstein study. 

5.1. Preservation of the Galactic Magnetic Field Direction 

Should the magnetic field as deduced from the Zeeman splitting of OH masers around massive 
SFRs show any Galactic-scale structure? This would require three important conditions to be 
satisfied. First, there must be a large-scale Galactic magnetic field, which is likely given pulsar 
(e.g., Manchester 1974) and extragalactic rotation measure observations (e.g., Sofue et al. 1986). 
Second, the OH clumps must "remember" the Galactic magnetic field after a magnetic compression 
of about three orders of magnitude from interstellar values. And third, intra-source reversals of 
the magnetic field in an individual massive SFR must be sufficiently few so as to not destroy a 
possible observed large-scale Galactic magnetic field by the inclusion of inferred local magnetic field 
directions that are opposite to the actual prevailing magnetic field direction near the source. 

Of these three conditions, the most difficult to justify on a physical basis is the assumption that 
the magnetic field direction at a site of massive star-formation would retain its initial (Galactic) 
direction through collapse. Ammonia observations of W3(OH) suggest that the density in OH 
maser clumps is n ~ 5 x 10 6 cm~ 3 (Reid, Myers, & Bieging 1987), which is a factor of 10 3-4 higher 
than found in giant molecular clouds and a factor of 10 6 greater than typical values in the ISM. It 
is unclear why the magnetic field in these clumps would retain the orientation of the field in the 
giant molecular cloud (GMC) during the compression phase, which may also entail many rotations 
of the clumps around the central condensation, and why the GMC would retain the orientation of 
the interstellar field. 
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5.2. Structure of the Galactic Magnetic Field 

Assuming that the second and third of the above conditions are satisfied, we can investigate 
whether there exists a large-scale Galactic magnetic field and what form it takes. The simplest 
imaginable Galactic field would be a circular field with no reversals in the portion of the disk where 
massive SFRs are found. Our 74-source sample discussed in §4.5 do not support this model. A 
plot of the predominant sense of the magnetic field direction (i.e., clockwise or counterclockwise 
as viewed from the North Galactic Pole) deduced from the OH masers near massive SFRs in 
this study is shown in Figure 3. For a circular magnetic field configuration in the Milky Way, 
the sense of the line-of-sight orientation (directed toward or away from the Sun) flips across the 
I = 0, 1 = 180° line. Noting this, 41 have magnetic fields directed in a clockwise sense as viewed 
from the North Galactic Pole and 33 have fields directed in a counterclockwise sense. This does 
not support a uniform circular field model without reversals. Indeed, it would not be inconsistent 
with a completely random distribution of magnetic field orientations. 

The rough equality of sources with clockwise and counterclockwise fields does not preclude an 
organized field structure. For instance, the magnetic field could be aligned along either concentric 
rings or spiral arms but with reversals between the rings or arms. Alternatively, the magnetic field 
as traced by interstellar OH masers could exhibit organized structure in small (kiloparsec-scale) 
patches, but not on a Galactic scale. Since our data do not imply a uniform sense of the Galactic 
magnetic field, a two-point spatial correlation function was calculated on the data in order to test 
the hypothesis that magnetic field directions of massive SFRs exhibit some degree of correlation 
over kiloparsec scales. We define the magnetic field orientation of two massive SFRs to be correlated 
if they are both oriented in a clockwise or counterclockwise Galactic sense. We similarly define a 
pair to be anti-correlated if the pair consists of one source with a clockwise field and one with 
a counterclockwise field. The correlation function was defined as the number of pairs of massive 
SFRs whose fields were correlated in the above sense, divided by the total number of pairs of 
sources. This function was calculated for a range of intersource distances. Since interpretation 
of the significance of the correlation function is difficult in isolation, it was compared against the 
same function evaluated for randomly-generated magnetic-field direction data at the locations of 
the massive SFRs in our sample. In the random sample, the positions of the maser sources were 
held fixed, but the line-of-sight field orientation of each source was randomly assigned as toward 
or away from the Sun. An ensemble of 10 5 such random field distributions was produced. The 
results are shown in Table 3. The percentages of trial runs showing a lesser and greater degree of 
correlation do not sum to 100% because some random distributions of magnetic fields produced 
identical correlation results. 

Since the data exhibit roughly equal likelihood of less and more correlation at the 10 or 20 kpc 
level, our data do not support the original hypothesis of Davies (1974) that there exists a uniform, 
clockwise Galactic field. On Galactic scales, pulsar rotation measure studies suggest that field 
reversals exist between spiral arms (e.g., Han et al. 2002; Rand & Lyne 1994). This would destroy 
any large-scale correlation on the Galactic scale, as spiral arms with oppositely-oriented magnetic 
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fields would be included at this scale. 

On sub-kiloparsec scales there are a number of source complexes in which two nearby massive 
SFRs have oppositely-oriented line-of-sight magnetic fields. The W49 complex provides a good 
example of this: G43. 148+0.015 (the Zeeman pair is between the sources W49A/R2 and R3 in the 
classification of De Pree, Mehringer, & Goss 1997) and W49S have fields directed primarily away 
from the Sun, while the field of W49N is directed toward the Sun. It is possible that this may 
result from a tangled field geometry associated with multiple cores in the same complex. 

Nevertheless, we do note a possible field organization on a ~ 2-kpc scale. Indeed, several 
patches of field coherence at this scale are noticeable in Figure 3. For instance, there is a cluster of 
sources with clockwise magnetic fields (i.e., oriented toward the Sun) near (X, Y) = (—4, —3) kpc, 
and another cluster with clockwise fields near (3, 5) kpc. 



The structure of the Galactic magnetic field can be studied with rotation measures of Galactic 
pulsars and extragalactic sources. We find that toward the second and third Galactic quadrants, 
which contain only six sources, all fields are oriented in the clockwise direction, while the nearest 
six sources in the first and fourth quadrants are consistent with a counterclockwise field. This 
agrees with rotation measure studies of the outer Galaxy, which also find an overall clockwise field 
(Lyne & Smith 1989; Brown & Taylor 2001). This also agrees with a reversal inward claimed in 
other papers from pulsar rotation measure estimates (Simard-Normandin & Kronberg 1979; Rand 
& Kulkarni 1989; Clegg et al. 1992). 

In order to compare our results with the Galactic magnetic field deduced from rotation mea- 
sures, we generated simulated rotation measure data from the line-of-sight magnetic field direction 
of our OH maser data. Our data consist of measurements of milligauss fields in situ around mas- 
sive SFRs, while extragalactic rotation measures are line-of-sight integrals predominantly through 
the microgauss fields of the interstellar medium. We calculated simulated rotation measures along 
rays of Galactic longitude within the Galactic plane using the line-of-sight field directions from OH 
masers. The rotation measure RM of radiation coming from a source at a distance D is given by 



where n e is the electron density in cm -3 , B is the magnetic field in gauss, and dl is measured in 
parsecs. Typical values for these variables in the interstellar medium are n e ~ 0.03 cm -3 (Lyne, 
Manchester, & Taylor 1985) and B\\ « 2 x 10~ 6 G (Manchester 1974; Simard-Normandin & 
Kronberg 1980; Heiles 1976). Given that magnetic fields implied by Zeeman splitting in massive 
SFRs exhibit some correlation on a scale of < 2 kpc (see §5.2), we postulated that the magnetic 
field around a massive SFR indicates the prevailing Galactic magnetic field direction in a sphere 
of radius r centered at the massive SFR. We adopted r ~ 0.5 kpc, which yields a contribution of 



5.3. Rotation Measure Comparison 



RM = 8.1 x 10 5 rad mT 2 




n e B • dl, 
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about 50 rad mT 2 to the rotation measure of a line of sight passing along the diameter of each 
sphere. Thus, we assigned a rotation measure of ±50 rad m -2 to each sphere, depending only on the 
line-of-sight magnetic field direction in the massive SFR, and summed the rotation measures of the 
spheres intersected by the ray. For simplicity, each sphere was assumed to contribute ±50 rad m~ 2 
to the total rotation measure, regardless of the path length of the ray within the sphere. Note that 
since the radiation propagates from D to 0, dl is oriented toward the Sun, and B • dl is negative 
when the magnetic field is oriented away from the Sun. Thus, a positive magnetic field generates 
a negative rotation measure. 

Our simulated rotation measures are plotted in Figure 4 for r = 0.5 kpc. Features of the 
plot are not highly sensitive to the value of r, but much lower values leave a large volume of 
the Galactic plane unsampled by our data, while much higher values result in a lot of overlap of 
the spheres. Superposed atop this plot are extragalactic rotation measures reproduced directly 
from Figure 9 of Rand & Kulkarni (1989), based on data from Simard-Normandin, Kronberg, & 
Button (1981). Our data agree reasonably well with the extragalactic rotation measure data. In 
the second and third Galactic quadrants, our data match the sign of the extragalactic rotation 
measures. While the extragalactic rotation measures are not generally of consistent sign in the 
first and fourth quadrants, agreement is found again at I = 15°, where the extragalactic rotation 
measures are decidedly positive. The simulated rotation measures we deduce from OH Zeeman 
splitting appear to fit the extragalactic data at least as well as, if not better than, the ring model 
of Rand & Kulkarni (1989). 



5.4. Spiral Arms 

In many models, one would expect the Galactic field to be organized along spiral arms. How- 
ever, investigating whether the Galactic field shows organization on a per-arm basis necessitates 
assignment of massive SFRs to spiral arms. Based on positions on both a plot of the Galaxy and 
on the longitude-velocity diagram in Figure 5, most of the sources in our study were assigned to 
a spiral arm, as listed in Table 4. The above correlation analysis was repeated on the subset of 
sources associated with each spiral arm, and the results are included in Table 3. For the most 
part, no obvious correlations exist except for sources assigned to the Norma arm, where 8 of 10 
sources (and all 7 in the fourth quadrant) have a magnetic field oriented clockwise with respect to 
the Galactic center. 

As noted above, it is difficult to make strong statements pertaining to field directions in spiral 
arms. The shape and locations of the spiral arms in the Milky Way are still a matter of debate, 
and many of the source distances used in this study are kinematic and therefore subject to large 
enough errors that some could be placed in the wrong arms. Even assuming that the near/far 
kinematic distance ambiguity has been resolved correctly for all sources, a typical distance error 
by this method is about 1 kpc. Furthermore, the OH maser velocity may differ from the velocity 
of the massive SFR, which may in turn deviate quite significantly from uniform circular motion, 
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as discussed in §4.1. Especially near I = 0°, where several spiral arms cross in a longitude- velocity- 
diagram, it is difficult to reliably assign massive SFRs to spiral arms. 

5.5. Intra-source Field Reversals 

The problem of relating OH Zeeman magnetic field measurements to the Galactic magnetic 
field can be inverted. Assuming that there exists a large-scale Galactic field, do our results allow 
us to comment either on whether the magnetic field threading OH maser clumps are related to the 
(pre-collapse) local Galactic field or on whether intra-source field reversals confuse the issue too 
much to see conclusive evidence of large-scale field organization? 

These are questions that cannot be answered conclusively with the VLA alone, due to its large 
beam size. For instance, no conclusive Zeeman pairs are found in the 1665 MHz transition in the 
VLA data for W3(OH). In contrast, Bloemhof, Reid, & Moran (1992) identified 23 Zeeman pairs in 
the same transition using VLBI. Significant blending of independent maser features within a single 
VLA beam often renders identification of Zeeman pairs ambiguous, as discussed in §4.5. When 
Zeeman pairs are identified with VLA resolution, they are often in short supply. We were able 
to find only one Zeeman pair in many of the massive SFRs shown in Table 1. One pair does not 
conclusively determine the prevailing magnetic field direction surrounding the region. Even in the 
simple case, depicted in Figure 6, of a conducting, spherical H II region in a uniform magnetic field, 
the line-of-sight projection of the magnetic field direction may vary depending on the inclination 
of the source (Bourke et al. 2001). A numerical simulation was run to estimate the fraction of 
sources whose line-of-sight magnetic field would be incorrectly inferred owing to projection effects 
for this simple geometry. Given an ensemble of ideal, spherical H II regions, each with a randomly- 
oriented prevailing magnetic field direction and a single magnetic field measurement taken locally at 
a random point on the projection of the sphere (as from a Zeeman pair), the line-of-sight direction 
of the measured field will be opposite to the line-of-sight direction of the prevailing magnetic field 
25% of the time. 

Intra-source field reversals have been seen across three of nine massive SFRs that were mapped 
in the main-line OH transitions with the VLBA: G351. 775-0.538 (Argon, Reid, k Menten 2002), 
W75N, and W75S (unpublished data). Thus, one-third of the population of massive SFRs mapped 
at milliarcsecond resolution show at least two distinct maser clumps with Zeeman pairs indicating 
opposite line-of-sight directions. While one might have expected that the magnetic field would be 
oriented randomly in a massive SFR owing to a complex evolutionary history that "scrambles" 
the pre-collapse field, massive SFRs with intra-source field reversals still exhibit an organized field 
structure. More likely, these reversals are a projection effect due to the geometry of the field lines 
relative to the observer. Still, the ionized sphere model is surely an oversimplification in some cases. 
For instance, the OH masers around W75N are aligned along two intersecting lines, one running 
roughly north-south and the other east-west (Baart et al. 1986). Our VLBA data show that the 
magnetic fields along the north-south line point consistently away from the Sun, while those in 
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the east-west line point consistently toward the Sun. Baart et al. suggest that the masers trace 
expanding shock waves, while Slysh et al. (2002) claim to have found evidence that the north-south 
structure is a thin, rotating disk. Whatever the explanation of this source may be, the observed 
intra-source field reversal is clearly not a simple projection effect in a static environment. 

5.6. The Random Component of the Galactic Magnetic Field 

Even if the Galactic magnetic field consists of a large-scale field with regular structure, the 
turbulent field component may hinder clear detection of the pattern. Synchrotron polarization 
observations suggest that the ratio of field strengths of the regular to turbulent magnetic fields may 
be between 0.6 and 1.0, with lower values more likely in spiral arms (Beck 2001). Assuming a 
uniformly random two-dimensional (in the Galactic plane) orientation of the turbulent component 
of the magnetic field, the line-of-sight orientation of the total field will be opposite to that of the 
regular component 8% to 20% of the time for the aforementioned range of field strength ratios. 
Thus, even if the pre-collapse magnetic field orientation is preserved at sites of OH maser emission, 
as many as one in five sources could imply a line-of-sight magnetic field opposite to the assumed 
underlying regular field. On the other hand, Brown &: Taylor (2001) present evidence that the 
random field may be partially aligned with the uniform field, rather than being isotropic. If this is 
indeed the case, the above estimates of the rate of discrepancy between the magnetic field direction 
implied by OH Zeeman splitting and the regular field direction may be overly pessimistic. Also, the 
random component of the Galactic magnetic field may become less significant between interstellar 
and GMC densities. 



6. Conclusions and Further Research 

We have found nearly 100 Zeeman pairs of ground-state OH masers using the data from Paper 
I. Combining these Zeeman pairs with others found in the literature, we have investigated the 
distribution in the Galactic plane of the line-of-sight orientation of magnetic fields in massive SFRs 
as deduced from OH maser Zeeman splitting. Our data do not support the original hypothesis of 
Davies (1974) that the large-scale magnetic field traced by OH masers is uniform in the clockwise 
sense, nor do we find that the inferred magnetic fields show a uniform direction in each spiral arm. 

All six sources in the second and third Galactic quadrants, which are typically 1 or 2 kpc 
distant, have a magnetic field oriented in the clockwise sense as viewed from the North Galactic 
Pole, while nearly all nearby sources in the first and fourth quadrants have a magnetic field oriented 
in the counterclockwise direction. This supports the hypothesis, deduced from pulsar rotation 
measures, that a reversal of the Galactic-scale magnetic field occurs near the galactocentric radius 
of the Sun. 

Our method for identifying Zeeman pairs in VLA maps of OH masers is limited by the blending 
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of maser components inevitable with a large (1" — 2") beamsize. The higher resolution of VLBI 
allows many more maser spots to be detected and Zeeman pairs to be identified unambiguously. 
The increased detail in VLBI images also allows for the identification of intra-source magnetic field 
reversals, the interpretation of which may be important both for understanding the conditions of 
the magnetic field in massive SFRs and for obtaining enough data to rigorously test the hypothesis 
that the magnetic field deduced from OH maser Zeeman splitting in massive SFRs is correlated 
with the Galactic magnetic field. To this end, we are in the process of imaging interstellar OH 
masers with the VLBA. The results will be published in a future paper. 



Table 2. Pairs with no Splitting 
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Av h 
(km s- 1 ) 


S b 

(Jy) 


fLSR b 
(km s- 1 ) 


Av h 
(km s- 1 ) 


A0 C 

(arcsec) 


G002. 143+0.010 




7.5 


1667 


2.09 


59.27 


0.34 


0.96 


59.27 


0.50 


0.026 


G009.622+0.195 




5.7 


1667 


4.67 


7.13 


0.29 


0.92 


7.14 


0.36 


0.032 


G043. 167+0.010 


W49 N 


11.5 


1665 


10.05 


7.88 


0.36 


32.08 


7.88 


0.43 


0.014 


G080.864+0.421 




4.1 


1665 


20.10 


-8.46 


0.26 


26.02 


-8.46 


0.24 


0.009 


G081. 721+0.571 


W75 S 


2.0 


1667 


0.74 


-1.12 


0.27 


1.56 


-1.14 


0.37 


0.070 


G081.871+0.781 


W75 N 


2.0 


1665 


26.00 


3.09 


0.35 


29.77 


3.07 


0.36 


0.002 










13.95 


0.63 


0.39 


17.31 


0.64 


0.39 


0.003 


Gill. 533+0.757 


NGC 7538 


2.8 


1665 


0.94 


-57.69 


0.30 


0.64 


-57.70 


0.35 


0.011 


G133. 946+1. 064 


W3 OH 


2.2 


1665 


32.94 


-47.45 


0.22 


20.61 


-47.45 


0.22 


0.015 


G350.011-l.341 




3.1 


1667 


2.25 


-23.74 


0.51 


2.60 


-23.78 


0.50 


0.016 


G351. 775-0.538 




2.7 


1665 


40.53 


-9.24 


0.25 


37.35 


-9.23 


0.27 


0.027 










11.66 


1.22 


0.50 


11.66 


1.22 


0.43 


0.099 



a Kincmatic distances assuming Rq = 8.5 kpc and ©o = 220 km s 1 , unless note in text says otherwise. 

b Peak flux density, center velocity, and line width (FHWM) of left- and right-circular polarized lines, respectively, as determined by 
Gaussian fit. 

c Angular offset between LCP and RCP lines. 
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Table 3. Correlation Results 



Separation 


Pairs of 


Percentage 


of Trial Runs Less/More 


Correlated 


(kpc) 


Sources 


Overall 


Car-Sgi 


: Cru-Sct 


Local 


Norma 


Perseus 


0.5 


43 


41/46 


14/65 


87/7 


.../... 


25/25 


0/63 


1.0 


111 


23/70 


56/31 


14/72 


13/38 


56/19 


0/63 


2.0 


324 


82/15 


81/15 


76/19 


44/25 


91/4 


19/44 


5.0 


1201 


52/45 


56/36 


66/28 


0/38 


95/2 


45/17 


10.0 


2197 


37/59 


42/48 


37/33 


0/38 


89/2 


30/45 


20.0 


2664 


63/36 


0/83 


37/33 


0/38 


89/2 


49/18 


Number of S 


ources 


74 


22 


17 


5 


10 


9 



Note. - - The figures represent the percentage of trial runs for which the correlation function 
evaluated on the random data were less and more correlated than the actual data. Results to any 
separation are cumulative; i.e., all pairs of sources with separation less than or equal to a given 
distance are included in the correlation function results to that distance. The total number of pairs 
of distinct sources with separation less than the figure in the first column is given in the second 
column. 
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Table 4. Sources and Field Directions by Arm 



Galactic 


Galactic 


Distance Magnetic 




Longitude 


Latitude 


(kpc) Ficld a 


References 


Carina-Sagittarius 


17.639 


+0.155 


2.1 


1 


20.081 


-0.135 


12.3 + 


1,7 


32.744 


-0.076 


2.3 


1,2 


34.258 


+0.153 


3.8 + 


1,2,3 


35.024 


+0.350 


3.1 + 


1,2 


35.197 


-0.743 


2.2 + 


1 


35.577 


-0.029 


10.5 


1 


43.796 


-0.127 


3.0 + 


1,2 


45.071 


+0.134 


4.7 + 


1,4,5 


45.122 


+0.133 


6.0 


1,2,4 


45.465 


+0.060 


6.0 + 


1 


49.469 


-0.370 


6.0 + 


1,6 


49.488 


-0.387 


5.5 + 


1,2,6 


285.26 


-0.05 


5.7 


2 


294.51 


-1.62 


1.4 


2 


311.64 


-0.38 


14.1 


2 


344.581 


-0.022 


2.0 


1 


345.505 


+0.347 


2.1 


1,2 


345.699 


-0.090 


1.6 


1,2 


348.549 


-0.978 


2.1 + 


1,2,7 


351.161 


+0.697 


2.3 


1 


351.416 


+0.646 


2.0 + 


1,2,7 


Crux-Scutum 


5.886 


-0.393 


3.8 + 


1 


11.03 


+0.06 


3.2 


2 


12.216 


-0.117 


3.3 


1 


12.680 


-0.181 


4.0 


1,8 


12.908 


-0.259 


4.0 + 


1 


28.199 


-0.048 


5.7 + 


1,3 


30.70 


-0.06 


5.5 + 


7 


323.459 


-0.079 


4.1 


2,10 


329.41 


-0.46 


4.3 


2 


331.34 


-0.34 


4.7 


7 


333.60 


-0.21 


3.5 + 


2 


339.62 


-0.12 


2.8 + 


2 


341.219 


-0.212 


3.2 


1 


343.128 


-0.063 


3.1 


1 


345.003 


-0.224 


2.9 


1,2,7 


350.011 


-1.341 


3.1 


1 


359.138 


+0.032 


3.1 + 


1 


Local 


69.540 


-0.976 


1.3 


1,3,9,11 


81.77 


+0.60 


2.0 


3 


106.80 


+5.31 


1.4 + 


3 


109.871 


+2.114 


0.7 + 


1,3,12 


213.706 


-12.60 


0.9 + 


1 


Norma 
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Table 4 — Continued 



Galactic 


Galactic 


Distance Magnetic 




Longitude 


Latitude 


(kpc) Field a 


References 




8.67 


-0.36 


4.8 + 


2 


9.622 


+0.195 


5.7 


1 


10.624 


-0.385 


4.8 


1 


330.95 


-0.18 


5.7 + 


2,7 


336.36 


-0.14 


10.3 + 


2 


336.83 


+0.02 


10.3 + 


2 


337.61 


-0.06 


12.8 + 


2 


337.71 


-0.05 


12.2 + 


2,7 


344.42 


+0.05 


13.4 + 


2 


350.11 


+0.09 


5.8 + 


2 


Perseus 


43.148 


+0.015 


11.5 + 


1,2 


43.165 


-0.028 


11.5 + 


1 


43.167 


+0.010 


11.5 


7 


70.293 


+1.601 


8.3 


1,3,11 


75.782 


+0.343 


5.0 + 


11 


80.87 


+0.42 


4.8 


3 


111.533 


+0.757 


2.8 + 


1 


133.946 


+1.064 


2.2 + 


1,3,7,13 


196.454 


-1.677 


2.8 + 


11 


Unassigned b 


0.547 


-0.852 


7.3 


1 


0.666 


-0.034 


8.2 + 


1 


0.670 


-0.058 


8.2 


1 


0.672 


-0.031 


8.2 + 


1 


40.42 


+0.70 


1.2 


1,2 


300.97 


+1.15 


4.4 + 


2,7 


305.81 


-0.24 


3.2 


7 


347.628 


+0.149 


9.8 + 


1,2 


354.73 


+0.29 


8.5 


2 


355.345 


+0.146 


23.1 + 


1 


359.436 


-0.103 


8.2 + 


1 



a Clockwisc (+) or counterclockwise (— ) magnetic field, as 
viewed from the North Galactic Pole. Note that "positive" mag- 
netic fields (pointing away from the Sun) are clockwise in the first 
and second Galactic quadrants but counterclockwise in the third 
and fourth quadrants. 

b Includes sources in the Galactic center. 

References. — (1) this paper; (2) Caswell & Vaile 1995; (3) 
Baudry et al. 1997; (4) Baart & Cohen 1985; (5) Zheng 1997; (6) 
Gaume & Mutel 1987; (7) Rcid & Silvcrstein 1990; (8) Zheng 
1991; (9) Dcsmurs & Baudry 1998; (10) Caswell & Reynolds 
2001; (11) in preparation; (12) Cohen, Brebner, k. Potter 1990; 
(13) Garci'a-Barrcto ct al. 1988 
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Fig. 1. — Sample spectra of windows used to find Zeeman pairs. Squares and triangles mark RCP 
and LCP emission, respectively. Gaussian fits to the lines are shown, and lines and features are 
numbered for reference. The number in the upper right corner marks the transition frequency, in 
MHz. (a-d) Windows in G351. 775-0.538, G032.744-0.076, K3-50, and G045.465+0.047. The 
1665 and 1667 MHz windows in (a) and (c) are centered at approximately the same position. The 
magnetic field and center velocity deduced from Zeeman splitting are shown in (a). 
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Fig. 2. — Histogram of the number of Zeeman pairs found with each magnetic field strength, 
binned by 1 mG. Positive (negative) field strengths represent magnetic fields oriented in the hemi- 
sphere away from (toward) the Sun. A best-fit Gaussian is also shown. The —25.4 mG pair in 
G081.871+0.781 is excluded for reasons outlined in §4.3. 
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Fig. 3. — Line-of-sight magnetic field directions deduced from OH maser Zeeman splitting. Seventy- 
four star-forming regions are plotted, distinguishing 41 with an overall magnetic field oriented in 
a clockwise sense from 33 with field oriented counterclockwise as viewed from above the Galactic 
center (i.e., North Galactic Pole). The locations of the Sun and Galactic center are marked. The 
spiral arms of Taylor & Cordes (1993) are indicated with shading and labelled by their common 
names. 



-35- 




Galactic Longitude (deg) 

Fig. 4. — Plot of simulated rotation measures generated from magnetic field directions deduced 
from OH Zeeman splitting. No simulated rotation measures are plotted for Galactic longitudes 
unsampled by the data of Figures 3 and 5. The dots represent extragalactic rotation measure data 
as taken from Rand & Kulkarni (1989) at Galactic latitude b = —15,0, and 15°. See §5.3 for 
further details. 
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Fig. 5. — Longitude- velocity diagram of magnetic field directions deduced from OH maser Zeeman 
splitting. Open and filled circles are defined as in Figure 3. 



Fig. 6. — Spherical bubble model of an H II region in a magnetic field. The field is oriented to 
the upper right and out of the page, inclined 30° to the line of sight. The arrows at a and b 
represent the local magnetic field at two points as might be deduced from Zeeman splitting. The 
line-of-sight projection of the magnetic field is out of the page outside the dotted curve (as at a) 
and into the page inside the curve (as at b). In the limiting cases of the magnetic field aligned 
along (perpendicular to) the line of sight, the line-of-sight direction of the field is reversed over none 
(half) of the projected area of the sphere. For a uniform distribution of magnetic field directions 
intersecting a sphere at all possible angles, on average the line-of-sight direction of the magnetic 
field is reversed over 25% of the area of the projected sphere. 



